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ABSTRACT 

Several photometric surveys for short-period transiting giant planets have targeted a number of open 
clusters, but no convincing detections have been made. Although each individual survey typically tar- 
geted an insufficient number of stars to expect a detection assuming the frequency of short-period 
giant planets found in surveys of field stars, we ask whether the lack of detections from the ensemble 
of open cluster surveys is inconsistent with expectations from the field planet population. We select a 
subset of existing transit surveys with well-defined selection criteria and quantified detection efficien- 
cies, and statistically combine their null results to show that the upper limit on the planet fraction is 
5.5% and 1.4% for 1.0 Rj and 1.5 Rj planets, respectively in the 3 < P < 5 day period range. For 
the period range of 1 < P < 3 days we find upper limits of 1.4% and 0.31% for 1.0 Rj and 1.5 
respectively. Comparing these results to the frequency of short-period giant planets around field stars 
in both radial velocity and transit surveys, we conclude that there is no evidence to suggest that open 
clusters support a fundamentally different planet population than field stars given the available data. 
Subject headings: planetary systems — open clusters and associations: general — techniques: photomet- 
ric 



1. INTRODUCTION 

The unexpected discovery of giant planets orbit- 
ing their host stars with periods of a few days 
([Mayor fc Quelodll995l iButler et alii 19971 ) made it clear 
that our theorie s of planet form ation and evolution are 
incomplete (e.g. lLin et al.l[l9 96). Although progress has 
been made toward understanding the the presence of 
short-period giant planets, many details of their forma- 
tion and evolution remain uncertain. 

That short-period planets form from circumstellar 
disks is clear, but the exact mechanism in which disk 
material condenses to form planets has yet to be es- 
tablished. Jupiter mass planets are thought to form 
through one of two mechanisms: core accretion (iMizunol 
119801 iBodenheimer fc PollacklU985 iPollack et al.lll996l) . 
in which a rocky or icy core reaches a critical mass and 
rapidl y accretes a n envelope of gas, or gravitational insta- 
bility (jBossI [T997I ) . in which the massive, cool disk frag- 
ments and collapses into giant planets. Core accretion 
theories have historically had difficulty producing planets 
of the correct masses and locations of the giant planets 
in our solar s ystem within a ~ 5 M yr protoplanetary disk 
lifetime fe.g. iPollack et afl 11996), although recent work 
has been more successful fsee iMovshovitz et al.l 120101 
and references therein). Gravitational instability, on the 
other hand, requires relatively cool, quiescent disk con- 
ditions and sufficiently short cooling timescales, both of 
which may not be physically real istic at the rad ii at which 
gas giants are thought to form (|Rafikovll2005l ). 

It is clear that short-period giant planets cannot form 
in situ, since protoplanetary disks at a few tenths of 
an AU are too hot and do not contain sufficient mass 
to produce a gas giant. Therefore, in addition to 
the uncertainties regarding the formation mechanism, 
we must face additional uncertainties related to the 
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means by which a planet migrates from its hypothe- 
sized birthsite in order to explain gas giants in short- 
period orbits. In what is labe l ed Ty pe I migration (fWardl 
119971 iGoldreich fc Tremaind [19801. the planet experi- 
ences torques from the massive gaseous disk, a loss of an- 
gular momentum, and subsequent inward motion toward 
the star. Sufficiently massive planets are capable of open- 
ing a gap in the circumstellar disk, which can also result 
in inward drift as the dis k accretes onto the star; this is 
called Type II migration (jLin fc Papaloizoulll986l: iWardl 
119971) . Recently, planet migration through planet-planet 
scattering and eccentricity pumping due to a stellar bi- 
nary companion coupled with tidal dissipation ("Kozai 
migration") have also been propos e d as migration mech- 
anisms (see [Nagasawa et all 120081 : IWu fc Murray! 120031: 
iMazeh et al.l 119971 ). The relative importance of each of 
these mechanisms to the process of planetary migration 
remains unclear. 

Each of these different formation mechanisms should 
imprint certain signatures onto the resulting planet pop- 
ulation. One hopes, with a sufficiently large observed 
sample of planets, to be able to determine the relative 
importance of each of these formation channels. From an 
observational point-of-view, the way to understand the 
origin of short-period planets is to search a large number 
of stars for planets and produce a correspondingly rich 
sample of exoplanets. The statistical properties of these 
planets and their observed frequency can then help to 
constrain models of planet formation and evolution. It is 
therefore essential that any sample of stars searched for 
planets has well characterized properties so that trends 
among the stars that are found to host planets can be 
r ecognized. 

Uanesl ([1996D suggested that open clusters were ideal 
targets for planet searches since they host such a uni- 
form and easily characterizable population of stars. In 
contrast to field stars, the metallicity, age, and distance 
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to cluster stars are relatively easy to determine. As such, 
open clusters represent a potential testbed for planet for- 
mation theories. The frequency of short-period planets 
provides information about how often gas giants form, 
migrate, and survive in a single, uniform, coeval pop- 
ulation of stars. Furthermore, observations of several 
clusters of different metalliciti es allow us to probe th e 
planet-metallicity correlation (|Fischer fc Valentil [2005). 
in which higher metallicity field stars are more likely to 
host a planet. 

Planet search surveys in clusters also allow us to char- 
acterize the degree to which the cluster environment it- 
self impacts the formation and survival of planetary sys- 
tems. The cluster properties at birth and throughout its 
evolution are critical to determining whether the cluster 
environment is capable of affecting the formation and 
survival of planetary systems. The degree to which the 
effects of close encounters and photoevaporating radia- 
tion present in dense stellar environments influences the 
formation, evolution, and survival of planets are not well 
understood. There is evidence that the observed stellar 
densities of old open clusters today are significantly less 
than they were at birth, and that the current environ- 
ments in old open clusters are not what they were during 
the first 5-10 Myr when giant planets were presumably 
forming. In particular, a given cluster may lose 10-80% 
of its stars in the process of emerging from its embedded 
stage, and thus only the initially most massive clusters 
survive and are observed as ol d open clusters today (see 
lLada fc Lad^l2tMlFrielll1^ . 

The degree to which planetary systems are af- 
fected by the cluster environment is still poorly un- 
derstood both observationally and theoretically. The- 
oretical i nvestigations of the importance of envi- 
ronment (iKobavashi fc Idal 120011: lAdams et~aTl 120061 : 
iBonnell et al.l 120011 : IArmitagel T2000'l suggest that open 
clusters are not dense enough to significantly affect the 
formation and survival rates of planets, especially within 
the ~ 5-30 AU in which gas giants that later migrate 
are thought to form. Ho wever, observations of the Orion 
Nebula Cluster (ONC) ([Eisner et al.1 120081 ) found that 
less than 10% of the stars harbor disks of mass compara- 
ble to that of the minimum mass solar nebula, and that 
the frequency of disks decreased with proximity to the 
massive central Trapezium stars. They also claim that 
the frequency of disks in the dense ONC is statistically 
different from that of the lower density Taurus cluster. A 
similar study of the Arches also detected a disk fraction 
that varies with p roximity to the massive, central stars 
(|Stolte et al.H2010D . 

Of the numerous planet detection methods avail- 
able to astronomers, photometric transit searches are 
among the most logica l choices for an open clus- 
ter survey (Ijanesl fl996l IPepper fc Gaudl [20051 [2001 
Ivon Braun et al.ll2005D . Transit surveys have the ability 
to monitor a large number of stars simultaneously, and 
benefit from the fact that cluster stars are concentrated 
over a small portion of the sky. Photometric surveys 
can also be performed on smaller, more easily accessible 
telescopes than compet i ng me thods of planet detection 
(jPepper fc Gaudil 120051 I2006D . The photometric data 
are simple to calibrate, and in the case of clusters the 
metallicities, distance, ages, and radii for all stars can be 
obtained with relative ease. Because the technological 



requirements for a successful photometric survey are ful- 
filled by modest and readily available instruments, tran- 
sit searches are a simple and direct method of detecting 
planets in these systems, with the caveat that the radial 
velocity follow-up needed to confi rm planet candidates i s 
not always as easily obtained (see lAigrain fc Pond feOO?). 

However, despite a number of transit s earches tar- 
getin g open clusters (|Hartman et al.l l2009t 
2006; M ocheiska et al.l 120081 1 20051 120061; 
2008: iBramich et al.l 120051: IBramich fc Hornd 1200. 
Hood et al.l 120051 : iBruntt et al.l 120031: IRosvick fe Rob 
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2006; Ivon Braun et all [2 005: 



Montalto et al.l l2007t iHowell et alj |2005l: IPepper et al l 



2008t Ivon Braun et all 12004 120051: | Lee et al.l 12004 
Hidas et al.l 120051: iStreet et al.l 120031: IRosvick fc Robbl 
120061 : iMontalto et al.1 120091 ) . no convincing planetary 
transits have been found. There are two possible 
expla nations for this paucity of planets ([Janes fc Kim 
2009): either there is something different about the 
planet populations of open clusters, or there are simply 
too few observed cluster stars to detect a planet from an 
otherwise "typical" planet popu lation. The frequenc y 
of short-period planets is ~ 1% ([Cumming et al.l [20081. 
and the geometric probability that such a planet transits 
is only ~ 10%, meaning that one naively expects to 
observe ~ 1000 stars before detecting even one transiting 
planet. Because the transit depths are small (~ 1%) 
and have a low duty cycle, both precise (to a few 
milimagnitudes) photometry and excellent temporal 
coverage are required in any transit search. When one 
considers the need for sufficient a signal-to-noise ratio 
and multiple observed transits for robust detection, it 
is generally the case that many more than 1000 stars 
must be monitored to detect any trans iting planets 
([Burke et all 120061 : Ivon Braun et al.l 120051 ). Since the 
richest open clusters contain a few thousand stars at 
best, the lack of detected transiting planets in any 
individual survey may simply be due to an insufficient 
number of target stars. 

While each survey alone may not observe a sufficient 
number of stars to have expected a planet detection, we 
estimate the total number of cluster stars observed in all 
transit surveys to be roughly 10,000, and that combined 
the sample of stars is rich enough to derive interesting 
upper limits on the frequency of short-period planets. By 
carefully statistically combining the null results of several 
photometric surveys, we provide a more stringent upper 
limit than is possible with individual studies, and then 
compare this limit to the frequencies of planets among 
field stars derived from both photometric and radial ve- 
locity surveys. As we will show, the typical detection 
efficiencies of these photometric surveys are insufficient, 
even with 10,000 total stars, to place tighter constraints 
on the frequency of short-period planets than surveys of 
field stars, and thus we conclude that these combined up- 
per limits are consistent with the observed frequency of 
short-period giant planets in the field. Given the avail- 
able data, we have no reason to suspect that the popula- 
tion of short-period planets in open clusters is anything 
other than "ordinary" . 

The paper is organized as follows: in |J5]we discuss the 
manner in which upper limits on planet frequencies are 
derived from transit surveys and how their individual re- 
sults can be combined, in fj3] we discuss our selection of 
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the transit surveys included in our analysis and the nor- 
malization of individual survey results. In ij5]we present 
our combined upper limits. $6] and $7] are devoted to 
discussion and our conclusions, respectively. 

2. METHODS 
2.1. Finding planets in a photometric survey 

The techniques used to find planets in any photomet- 
ric survey are broadly similar. All surveys attempt to 
achieve high photometric accuracy with maximal tem- 
poral coverage over the longest possible period of time. 
The better the temporal coverage, the more sensitive a 
survey is to all transits, and the longer the duration of 
the survey, the larger the range of plane tary periods that 
are de tecta ble. We refer the reade r to iPepper fc Gaudil 
(pOOfih and Ivon Braun et al.1 (|2005| ) for extensive discus- 
sions of the factors that determine a survey's sensitivity 
to transits and the design of successful surveys. 

Once the data have been collected, one produces 
light curves and searches for periodic variability. The 
box-fitting least-squ ares (BLS) method introduced by 
iKovacs et al.l <|2002f ) is a popular means to do so within 
the transit community. The algorithm searches for peri- 
odic, rectangular deviations from a flat light curve, and is 
an objective, repeatable means of identifying transit can- 
didates. Conversely, some authors may choose to identify 
transits by eye. Regardless of the method used to identify 
photometric eclipses, the planetary nature of convincing 
transit candidates must then be confirmed with radial 
velocity follow-up observations, since astrophysical false 
positives can mimic the transit of a Jupiter-sized object. 

In the case of a null result, some authors choose to 
carefully quantify their detection efficiencies and place 
upper limits on the frequency of short-period planets. 
One typically injects a large number of simulated, limb- 
darkened transits into constant simulated stars or light 
curves from the survey itself, and then subjects these 
transits to the same detection algorithms with the same 
selection criteria as the original data. The fraction of in- 
jected transits that are recovered quantifies the detection 
efficiency. This efficiency is a complicated combination of 
effects that is discussed in more detail in the next section. 

While the general process is similar for all surveys, the 
details of how each author chooses to perform the steps 
varies slightly, and in practice it is not necessarily cor- 
rect to directly compare or combine the results of two 
surveys. Authors may make different assumptions about 
the period distribution of planets, or quote their results 
for differing planetary radii, etc., all of which are im- 
portant when comparing the results of multiple different 
surveys. In order to combine the results of several sur- 
veys, we must re-normalize, to the extent that we can, 
to a common set of assumptions. 

2.2. Quantifying upper limits 

In this section we discuss the mathematical description 
of the detection efficiency of a transit survey, the calcu- 
lation of upper limits on the frequency of planets using 
a null result, and the method for combining several nor- 
malized surveys into a single upper limit. 

In general, the expected number of planets detected in 
the radius range R p to R v + dR v and peri od range P to 
P + dP can be written (jBurke et al.ll2006j) 



d 2 N e 



dRpdP 



d 2 P 
dR p dV 



N, 



kVtr,kVdet,k, (1) 



where k is the star, the sum is over all N* stars, and 
r Pmem,k is the probability that the star is a cluster 
member. f p is the fraction of stars that host short- 
period planets distributed as d 2 P/dR p dP. We assume 
d 2 V I dRpdP is independent of the stellar mass. This dis- 
tribution is poorly constrained by exoplanet statistics, 
and is usually assumed as a prior. We will assume that 
this distribution is a delta function in the radius and 
uniform in the logarithm of the period: 
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where R p is the planetary radius of interest. Vt r ,k is the 
geometric probability that a planet transits its host star, 
such that 



Vtr,k — 



Rp + R* 



txM* 1/3 (R p + R+)P 
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where i?* is the radius of the star and a is the semima- 
jor axis of the orbit, where we have assumed the orbit is 
circular. Vdet,k is the probability that a planet of radius 
Rp and period P will be detected around star k if it tran- 
sits, given the precise temporal coverage, precision, and 
signal-to-noise ratio of the observations. In practice Vtr 
is trivial to calculate, while Vdet is emphatically not so. 
It is, in general, the single most difficult factor to char- 
acterize in the entire formalism because it represents a 
complex interplay between numerous different observa- 
tional effects. Simple analytic estimates of Vdet tend to 
overestimate the probabilit y that the survey will be ca- 
pable of detecting transits ()Beattv fc Ga udi 2008) . The 
best way to characterize Vdet is to perform Monte Carlo 
simulations a posteriori with transits injected into the 
light curves of constant stars observed in the survey or 
simulated constant stars. 

When authors choose to perform an analysis of the 
detection efficiencies with Monte Carlo simulations, the 
resulting efficiencies generally account for Vdet, Vtr, and 
d 2 V I dRpdP . By injecting transits with randomly sam- 
pled parameters from an assumed period distribution 
into stars in the survey that were searched for transits, 
one automatically accounts for the intrinsic mass func- 
tion of stars and planet period distributions. Factors 
such as noise present in the light curves and the ob- 
serving window, which are crucial to determining Vdet, 
are also automatically included in transit injection and 
recovery schemes. V mem is generally determined sepa- 
rately through proper motion measurements, proximity 
of a star to the cluster main sequence, off-cluster com- 
parison fields, or modeling of the field star population 
to determine contamination levels. After the detection 
efficiencies are quantified, one has 
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where N exp is the number of transiting planets that are 
expected to be detected, and f p is the only unknown fac- 
tor. By assuming that Vdet,k is independent from Vmem.k 
for each star we can define an average membership prob- 
ability (Vmein)- While these two factors are not strictly 
uncorrelated, it is a reasonable approximation and sim- 
plifies the analysis considerably. Using this assumption 
we can average over all of the stars in a survey such that 

N exp = NJ p {V mem ){V e ), (5) 

where (V e ) is the average probability over all stars that 
a planet will both transit and be detected. 

In the case in which no planets are detected over the 
course of the survey, one can derive an upper limit on the 
frequency of planets of a given radius within a specific 
period range. Formally, the probability of seeing Ndet 
successful detections in ./V* independent trials, each with 
a probability of success f p {V e )(V m em) is given by the 
binomial distribution 



V(N det ;N^f p (V e )(V mem )) 

' N* ' 



(6) 
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Because the number of stars (and therefore trials) is large 
and the probability of detecting a planet around a given 
star is small, we can assume the Poisson approximation 
to the binomial distribution without compromising our 
results. The probability of detecting N de t events when 
the number of expected events is N exp is 



]\fN iete -N e 

P(N det ;N exp )= 



N det \ 

In the case of a non-detection, Ndet — and 



V(0; N exp ) = 



(7) 



(8) 

As mentioned earlier, the number of detections we expect 
for a given survey goes as 

N exp = f p N+ (T mem )(V e ) = f P N cl {V t ) , (9) 

where N c i = AT* (V m em ) is the number of cluster members 
observed. The 95% confidence upper limit on the planet 
fraction (when V(0; N exp ) — 0.05) is then 



fp < 
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The multiple surveys that we utilize in this paper each 
quote a value of f p or N c i(V e ) specific to their obser- 
vations, chosen planetary radii, and period ranges. Pro- 
vided the limits or detection efficiencies for all surveys are 
quoted for the same planetary radii and period ranges, 
the upper limits from individual surveys can be combined 
using 

-In (0.05) 



< 



where f to t is the combined upper limit on the planet frac- 
tion, and s indexes individual surveys. The expression 
for ftot can be written simply as 



ftot < 



(12) 



3. SELECTION OF TRANSIT SURVEYS AND BASIC 
ASSUMPTIONS 

Although the general planet search techniques and sub- 
sequent quantification of detection efficiencies are rel- 
atively similar among the photometric surveys consid- 
ered here, there are always slightly different assumptions 
made during the process by different authors. In the 
transit surveys we ultimately selected for our study, the 
two main areas in which differences were apparent were 
the assumptions made about the period distribution of 
planets and the estimates of cluster membership. Such 
differences are relatively straightforward to re-normalize, 
and do not pose a fundamental obstacle to combining the 
surveys. However, in the cases in which authors do not 
quantify their detection efficiency sufficiently, or do so 
only in terms of the most simplistic analytic estimates, 
we cannot re-normalize their results, and therefore can- 
not include them here. We also have no way to quantify 
the detection efficiency of transit surveys in which the 
transits were identified "by eye", and so these too can- 
not be utilized. 

We choose transit surveys with the following character- 
istics that make it possible to re-normalize and combine 
their results: 

1. Transits are detected through the use of an auto- 
matic transit detection algorithm with rigorous, de- 
fined thresholds of detection. Specifically, "by eye" 
transit detection is not sufficient, since detection 
probabilities are not quantifiable. 

2. A simulation of detection efficiency was conducted 
using constant stars, and injected transits are limb- 
darkened and have a representative sampling of pe- 
riods, inclinations, and planetary radii. This allows 
for a realistic quantification of elements such as the 
window of observations, instrument noise, and the 
effectiveness of the transit detection algorithms. 

3. Survey results must be accompanied by an estimate 
of cluster membership probabilities. It is essential 
that we have a reasonable estimate of the actual 
number of observed cluster stars, as opposed to a 
blind mixture of cluster and field stars. Without 
a membership estimate, any upper limits derived 
from the survey results for the frequency of planets 
in clusters is not meaningful. Surveys are included 
in cases in which membership is not explicitly esti- 
mated but sufficient information to do so is avail- 
able. 

Our analysis therefore includes the results from 6 
separate transit surv eys of 6 op en clusters: M3 7 
(IHartman et al.|[2009l). NGC 124 5 (|Burke et al.|[200l , 
NGC 188 (IMocheiska et"al|[200l hereafter M08), NGC 
6791 (IMocheiska et al.ll2005L hereafter M05), NGC 2158 



(11) (IMocheiska et al 



20061 hereafter M06), and NGC 2362 



(jMiller et al.ll2008l) . In the case of NGC 6791, three sepa- 



rate surveys targeted the cluster and we choose to utilize 
only one for simplicity. The fundamental cluster param- 
eters (distance, age, etc.) of the systems observed in the 
relevant surveys are given in Table [T] and a discussion of 
the individual surveys follows in g] 

In addition, we adopt the following conventions to nor- 
malize the results from the six surveys: 
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TABLE 1 

Fundamental cluster parameters 



Cluster (a°,<5°) (l,b) distance (pc) age (Myr) [Fe/H] References 3 - 

M37 88.074 +32.553 177.637, +3.094 1500 + 100 550 + 30 0.05 + 0.04 2,2,2 

NGC 188 12.108 +85.255 122.865, +22.384 1710 + 80 6000+ 1000 -0.04 + 0.05 7,8,6 

NGC 1245 48.68 +47.25 146.64, -8.92 2800 + 200 1040 + 90 -0.05 + 0.09 1,1,1 

NGC 2158 91.054 +24.097 186.634, +1.781 3600 + 400 2000 + 300 -0.06 5,5,10 

NGC 2362 109.65 -24.98 238.20, -5.58 1480 5±^ •■• 3,3,- 

NGC 6791 290.221 +37.772 69.958, +10.904 4200 8000 + 500 0.39 + 0.01 4,4,9 



a R eferences for cluste r par ameters, given in the order: distance, age, mctallicity. Refere nces are as fol lows: 
1) IBurke et al.l J2004T) . 2.) IHartman et al.l j200Sl ) , 3.) iMoitinho et al.l ||200l) ) 4.) IChabover et al.l 1T999T), 5. ) 
ICarraro et al.l [120021) .6 ) Ivon Hippel & Saraic dini (19981). 7.1 ISaraiedini et alj l[l9991) . 8.) IDinescu et al.l H1995). 
9.) ICarraro et al.l 12006T) . 10.) IJacobson et al.1 H2"00Sh 



1 . We assume that the underlying planet population is 
an even logarithmic distribution in period. In cases 
where the authors of the papers have not made a 
similar assumption, we correct for the difference 
between the assumed distributions. 

2. We quote upper limits at the 95% confidence level, 
unless otherwise stated. 

3. "Hot Jupiters" (HJs) are defined to be planets with 
periods in the range 3 < P < 5 days, and "Very 
Hot Jupiters" (VHJs) those that have periods of 
1 < P < 3 days. 

3.1. Planet Candidates 

Some of the selected surveys detected transit-like 
events of the appropriate depth among their light-curves, 
which represent candidate transiting planets. We assume 
here that all candidates are false positives. We justify 
this assumption based on the fact that most candidates 
are revealed as false positives upon c loser inspection (see 
IBrowrJ [20031 : lEvans fc Sack^l2lTI(it ). none of these can- 
didates have been confirmed as planets through radial 
velocity follow-up observations, and, as we will describe 
later in our upper limits are consistent with radial 
velocity and transit surveys of field stars even without 
the inclusion of the candidate planets (and thus our as- 
sumptions are conservative). 

4. INDIVIDUAL CLUSTER SURVEYS 

The upper individual limits derived for each cluster 
and the details of the re-normalization of each survey 
follow. Our final adopted upper limits derived from in- 
dividual surveys are listed in Tabled] 

4.0.1. MS7 

The IHartman et~aTl (2009) survey of M37 utilized in- 
jected, limb-darkened transits with randomly selected 
periods, radii, and inclination angles to arrive at their 
detection probabilities, and determined cluster member- 
ship by using the luminosity function of a nearby blank 
field and a narrow strip on the CMD enclosing the clus- 
ter main sequence to arrive at membership probabilities 
as a function of r magnitude. Planets are assumed to 
be evenly distributed in logP and cosi. We adopt the 
authors' quoted upper limits without adjustment. 

4.0.2. NGC 188, NGC 2158, and NGC 6791 



The M08 survey of NGC 188, M06 survey of NGC 
2158, and M05 survey of NGC 6791 all utilize identical 
transit injection and recovery schemes to ascertain their 
detection efficiencies. A total of 432,000 limb-darkened 
transits are injected into observed light curves without 
transits for each survey and run through the detection 
algorithms for recovery. They inject signals correspond- 
ing to periods from 1.05 to 9.85 days in linear steps of 0.2 
days, radii from 0.95 to 1.50 Rj in steps of 0.05 Rj, and 
cos i from 0.0125 to 0.9875 in steps of 0.025. Detections 
are defined as "marginal" or "firm" by the authors if they 
fulfill one or both of the transit selection criteria: the pe- 
riod of the injected transit was recovered to within 2% of 
the true value or that of an alias, and certain thresholds 
were met for the values of the BLS statistics. We con- 
sider both classes of detections here, but use the "firm" 
detections in all combined upper limits quoted later in 
the text and in all plots, since a "marginal" detection can 
represent an event with only a single detected transit. 

There are two ways in which we must adjust the pub- 
lished values in order to compare them to other surveys 
in our sample. The authors assume that planets are dis- 
tributed linearly in period, while we have adopted the 
convention that the distribution is uniform in the log- 
arithm of the period, and so we must rescale their re- 
sults to reflect our choice. Secondly, estimates of cluster 
membership are lacking or in need of fine tuning, and we 
therefore estimate membership probabilities based on the 
information provided in each of the three papers. 

In all three papers the authors assume that planets are 
uniformly distributed in period. The average probability 
of detecting a planet of any of the considered radii (0.95 
to 1.50 Rj) and inclinations (cosi from 0.0125 to 0.9875) 
is given by the period-frequency histograms in Fig. 4 
in M08, Fig. 6 in M06, and Fig. 5 in M05. In order 
to rescale the detection efficiencies to a planet distribu- 
tion uniform in the logarithm of the period, we apply 
a multiplicative scaling factor to each period-frequency 
histogram bin. This scaling factor, Sp dist , is given by 

P — P 

1 max 1 min /1 o\ 

Z>Pdi3t,j - T7, \~' ^ l6 > 

Pbinj In ( jP^ J 

where j is the bin index, P max (P m in) is the largest 
(smallest) period in the period range of interest (i.e. 1 
and 3 days, corresponding to HJs with 1 < P < 3 days), 
and Pbin,j the period upon which the histogram bin is 
centered. The net effect of this correction is to increase 
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TABLE 2 

Upper limits on the frequency of HJs and VHJs 











VHJ 






HJ 




Cluster 


N cl 




1.0 Rj 


1.2 Rj 


1.5 Rj 


1.0 Rj 


1.2 Rj 


1.5 Rj 


M37 


1450 




0.027 


0.025 


0.023 


0.083 


0.077 


0.069 


NGC 188 


797 


firm 


0.222 


0.160 


0.067 


1.87 


1.34 


0.560 






marginal 


0.051 


0.042 


0.027 


0.178 


0.144 


0.093 


NGC 1245 


870 




0.240 


0.170 


0.064 






0.360 


NGC 2158 


2460 


firm 


0.076 


0.049 


0.009 


0.394 


0.255 


0.047 






marginal 


0.021 


0.015 


0.006 


0.052 


0.038 


0.016 


NGC 2362 


475 




0.521 


0.369 


0.142 


4.99 


3.30 


0.775 


NGC 6791 


1997 


firm 


0.082 


0.053 


0.008 


0.346 


0.221 


0.033 






marginal 


0.025 


0.018 


0.006 


0.069 


0.049 


0.018 



■ Please note that AT C ; refers to the estimated number of cluster members observed with 
sufficiently precise photometry to detect planets. 



the detection probability of short-period planets relative 
to that of longer period planets. 

Each period-frequency histogram bin is an average over 
all planetary radii and inclinations. If we wish to calcu- 
late the probability of detection for a particular radius, 
we must multiply the bins by yet another scaling factor 
that quantifies whether a transit of a particular plane- 
tary radius is detected more or less efficiently than av- 
erage. For this adjustment we use the radius-frequency 
histograms in Fig. 5 M08, Fig. 8 in M06, and Fig. 7M05, 
which plot the average probability of detecting a planet 
of a particular radius, where the probabilities represent 
an average over all considered periods and inclinations. 
The scaling factor is given by 



(V)' 



(14) 



where Pr v is the probability, read from the radius- 
frequency histogram, of detecting a planet of radius R p , 
and (V) is the average detection probability over all radii. 

With these scaling factors we can then calculate the de- 
tection efficiency for any desired period range and planet 
radius, T> e from Eqn. [5j described as 



N hi ,. 



P, 



max ± tmn 



P-n 



(15) 



where Vj is the probability of detection in each period bin 
in the frequency-period histogram, Nu ns is the number 
of bins contained in the period range P m in < P < P ma x, 
and AP is the period width of each individual bin. This 
technique is applied to all three of the clusters NGC 188, 
NGC 2158 and NGC 6791. 

4.0.3. NGC 6791 membership 

As noted earlier, it is critical that we have estimates 
of the actual number of cluster stars observed. The M05 
paper provides estimates the number of cluster and field 
stars in the sample used for the transit analysis. Candi- 
date cluster members are selected as those stars that lie 
within 0.06 mag of the V — Rvs.R cluster main sequence. 
While this significantly reduces the field star contamina- 
tion, there is still a non-negligible number of field stars 
among these candidate cluster members that happen to 
have the same colors as the cluster main sequence. To 



quantify this contamination we use t he BesancQn stellar 
population synthesis galaxy model (R obin et al.1 12003ft 
to estimate the number of field stars that should be 
present in the region of the cluster CMD used to de- 
termine membership. Secondly, to confirm the result of 
the first method, we then estimate the typical field star 
density on the CMD and check that the resulting mem- 
bership is consistent with the estimates derived from the 
galaxy model. 

We first utilized the Besancgn galaxy model to esti- 
mate the number of field stars expected in the M05 field 
of NGC 6791. The models were run for a field of view of 
0.144 deg 2 with 25 dust clouds evenly spaced every 2pc 
from the observer, each with an extinction o f 0.019 for a 
total Ay = 0.477 from iSchlegel et al.1 (|1998ft . We found 
that this yielded more realistic star counts than using a 
uniform diffuse extinction over the entire distance to the 
cluster. We restrict the sample in the same manner as 
M05 to stars with R magnitudes 17 < R < 19.8, which 
corresponds to all stars below the main sequence turnoff 
and an rms uncertainty in the R magnitude in the M05 
survey of less than 0.05 mag. All stars from the Besancgn 
simulation within 0.06 mag of the cluster main sequence 
were then selected as field star contaminants. 1523 of 
the 5225 field stars predicted by the model were within 
this region. The M05 member selection method therefore 
would have found 1523 cluster and 3702 fields stars in this 
model data. It is important to note that we do not trust 
the absolute star counts produced by the model, but ex- 
pect the fractional number of stars that have colors simi- 
lar to the cluster main sequence to be more accurate. The 
ratio of stars detected off of the main sequence in M05 
and the simulation is 2871/3702 « 0.78, and thus M05 
sees ~ 22% fewer stars that are more than 0.06 mag from 
the MS than the simulations of the same field. We scale 
the number of stars present in the model near the MS 
to the ratio of the number of stars detected in the field 
portions of both the simulated and actual cluster CMD. 
We would therefore expect that 1523 x 0.78 = 1181 stars 
selected as cluster members in M05 are actually field star 
contaminants. This implies that only 1997 of 3178 can- 
didate cluster stars are likely to be true members, so 
N cl = 1997. _ 

Another simple, and somewhat crude, method of esti- 
mating the field star contamination is to simply find the 
average density of stars on the cluster CMD near, but 
not on, the main sequence. This density, multiplied by 
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the area covered by the main sequence (in mag 2 ), yields 
an estimate of the number of field stars contaminating 
the selected sample of cluster stars. We find surface den- 
sities of 650-700 stars mag -2 and thus the main sequence 
area of 0.12 x 3.54 mag 2 yields roughly 300 contaminat- 
ing stars along the MS strip on the single-field CMD 
in Figure 2 of M05. This would correspond to roughly 
1200 stars over all four observed fields, and about 2000 
true cluster members. While the method is approximate 
at best, it yields a reassu ringly similar result. Finally, 
iKaluznv fc Udalskil (|1992f ) find that the cluster CMD is 
subject to 30% contamination by field stars for stars with 
V < 20 mag in NGC 6791, which is again comparable to 
our estimates. 

4.0.4. NCG 188 membership 

M08 does not provide any estimate of the num- 
ber of cluste r stars observed in NGC 188. However, 
iPlatais et al.l ()2003l . hereafter P03) provides a catalog 
of stars down to V = 21 in NGC 188, and obtained 
proper motions and membership probabilities using mea- 
surements from earlier epochs taken with photographic 
plates. These data cover a magnitude range in V simi- 
lar to that of the stars utilized for the transit search in 
M08. Since the P03 field only covers a fraction of the 
field of the M08 photometric search, we use the central 
field (labeled A3 and Bl in M08) of the M08 survey to 
derive estimates of the field star density across all 6 of 
the other observed fields. This reference field is centered 
on the cluster and ~ 450 stars are detected in M08. This 
same field in the P03 survey contains 612 stars, all but 
6 of which have measured individual cluster membership 
probabilities. Of these ~ 600 stars, 372 stars have mem- 
bership probabilities P mem > 50%, and 234 stars have 
Pmem < 50%. To be conservative, we assume that all 
stars with P mem < 50% are field stars. For a reference 
field size of 11'. 4 x 11'. 4 this implies a field star density of 
~ 6500 field stars deg~ 2 down to V ~ 21. However, since 
M08 detected fewer stars in this same field, we scale the 
P03 field star density by the ratio of stars detected 
in the M08 and P03 surveys respectively. This gives a 
field star density of |f| X 6482 stars/deg" 2 = 4883 field 
stars deg -2 . M08 covers 7 11'. 4 x 11'. 4 fields, so we ex- 
pect 1234 field stars. Since 2031 stars are detected on 
all 7 chips combined, this suggests that 1234 are field 
stars, and the remaining 797 are cluster members, and 
so N c i — 797. We note that the earlier estimates that 
made use of the Besancgn model are not feasible here, 
since we have no estimate of the number of field stars 
with which to calibrate the star counts returned by the 
model. 

4.0.5. NGC 2158 membership 

The M06 analysis of NGC 2158 also did not include any 
estimate of the number of cluster members observed. In 
this case, existing surveys of the cluster contained either 
no membership information, covered too small a field 
of view to be easily comparable, or did not cover the 
relevant magnitude ranges. To be conservative, then, 
we make the assumption that M06 field 1 centered on 
a w 91.65°, S R3 23.9° contains only field stars. This is a 
11 .4 x 11 4 region centered roug hly 16' from the clus- 
ter center. iKharchenko et ahl ([1997D estimates the max- 



imum radius of the cluster to be ~ 15', so this field 
likely contains at least some cluster stars, although the 
st eep radial density profil e of the cluster (see Fig. 6 
in IKharchenko et all 119971 ) suggests it will be minimal. 
Note also that field 1 in M06 is also the most distant field 
of the four from nearby M35, which produces a slight gra- 
dient in the star counts across the M06 fields. M06 finds 
675 stars in this field with photometry of suitable quality 
to be used in the search for transits. If we assume all of 
these stars are field stars it implies that there are a to- 
tal of 2700 field stars across all four fields. A total 5159 
stars, both cluster and field, with suitable photometry 
were found in all four fields. We can then estimate N c i 
as 5159 — 2700 ~ 2460. This is a conservative estimate 
of the fraction of the observed stars that are true cluster 
members, since it is unlikely that all 675 stars were truly 
field stars. This approach preserves the fidelity of the 
upper limits we derive using this value. 

4.0.6. NGC 1245 

The iBurke et~aTl pOOl survey of NGC 1245 includes 
detailed detection probability simulations and careful 
treatment of cluster membership. Membership proba- 
bilities are determined using a star's proximity on the 
color magnitude diagram of NGC 1245 to the best fit 
isochrone. Because each individual star is assigned a sep- 
arate membership probability, we do not correct further 
for field star contamination as we did in the M05 NGC 
6791 survey. Planets are assumed to be uniform in the 
logarithm of the semimajor axis, which is equivalent, up 
to an irrelevant normalization constant, to a distribution 
uniform in log P. As such, the upper lim its on the planet 
fraction as derived by IBurke et al.l ()2006f ) are used in our 
analysis without adjustment, save for upper limit on 1.5 
Rj companions in the HJ period range, which we discuss 
b elow. 

Bur ke et"aT1 (|2006l ) takes HJs to the be planets with 
periods 3.0 < P < 9.0 days, while we have assumed the 
the HJ period range to be 3.0 < P < 5.0 days. The 
authors quote an upper limit f p < 36% for the period 
range 3.0 < P < 6.0 days for 1.5 Rj planets. One can 
see from the middle panel in Fig. 8 of their paper that 
the upper limits for a period range of 3.0 < P < 5.0 days 
will be negligibly different than the quoted values for 
3.0 < P < 6.0 days. If anything, the constraints on the 
slightly longer period range should be weaker than those 
of the 3 < P < 5 day range, and so we are conservative in 
adopting the value of f p < 0.36 for the 3 < P < 5 day 1.5 
Rj planets. No limit could be placed on the frequency 
of 1.0 Rj planets in this period range. 

4.0.7. NGC 2362 

As with the other surveys, the IMiller et ahl (|2008f ) 
analysis of NGC 2362 light curves relies on objective 
and clearly defined thresholds of detection and a tran- 
sit injection scheme to characterize the sensitivity of the 
survey. A r ough estimate of cluster membership from 
llrwin et al.l (|200 8i) using both field star counts from the 
Besancgn galaxy model and simple radial distribution ar- 
guments arrives at a field star contamination of roughly 
60% over the magnitude range of stars utilized in the 
transit search. 

The authors quote upper limits on planet frequency for 
the 1 < P < 3 days and 3 < P < 10 days period ranges, 
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while here we are interested in the 3 < P < 5 range. We 
use their Figure 7 to estimate the fraction of recoveries 
in the 3-10 day range that would also be detected in the 
3-5 day range. Using both the 2 transit and 3 transit 
curves we estimate that / 3 _ 5 k 0.58 x /3_io- Combined 
with the number of expected detections in their Table 
7, this corresponds to f p < 4.99 for 1.0 Rj planets and' 
fp < 0.76 for 1.5 Rj planets at 95% confidence in the 
3-5 days period range. Note that f p > 1.0 simply means 
that the data is only able to place limits on the fraction of 
stars with more than one planet in a given period range. 

5. UPPER LIMITS ON THE PLANET FREQUENCY 

With the upper limits on the planet frequencies for 
each of the transit surveys in hand we can provide a 
combined upper limit using Eqn. [T^Jin Upper limits 
are presented for planetary radii of 1.0 and 1.5 Rj in the 
H J and VH J period ranges here and in Table |3J We find 
upper limits to the fraction of planets with radii of 1.0 
Rj and 1.5 Rj to be 5.5% and 1.4%, respectively. We - 
find upper limits on the frequency of VHJ of 1.4% and 
0.31% for 1.0 Rj and 1.5 Rj respectively. 

In general, the average radius of short-period plan- 
ets is greater than 1.0i?j and less than 1.5Rj, with 
mean and median in the range ~ 1.2 — 1.3 Rj (from 
www.exoplanets.org, as of June 2010), and thus to com- 
pare to field surveys for planets, we would prefer to 
quote upper limits for this radius range explicitly. How- 
ever, very few authors explicitly calculate upper limits for 
these radii, and thus we must interpolate. If we assume 
that the detection probability increases linearly with in- 
creasing planet radius (as Figures 7, 8, and 5 in M05, 
M06, and M08, respectively imply), we can also estimate 
upper limits for 1.2 Rj objects. For this linear interpo- 
lation, the upper limits are 0.95% and 3.9% for the VHJ 
and HJ period ranges, respectively. 

Unfortunately, the vast majority of surveys fail to pro- 
vide errorbars on their upper limits or expected number 
of detections. In addition, the nature of our membership 
probabilities are such that they represent estimates only. 
For both of these reasons we do not quote uncertainties 
on our upper limits. 

5.1. Metallicity considerations 

There is a well known correlation between the metal- 
licity of the host star and the likelihood that the star 
bears a planet that can further inform our analysis. 
IGonzalea (|1997[ ) first noted that planets tend to be ob- 
served around stars with higher than average metallic- 
ity. If this correlation is due to the fact that primordial 
clouds of higher metallicity more efficiently form planets, 
rather than p ollution due to the proc e ss of planet for - 
mation itself (jSantos et alJl2000T, [20011: lLaughlin|[2000l ). 
we would expect to see relatively more planets in open 
clusters of higher metallicity, and give more weight to the 
null results of transit sur veys of more metal-rich clusters. 
iFischer fc ValentH ([2005) provides the em pirical relation- 
ship (for a recently updated version see Uohnson et al.1 
120101) between the probability of hosting a planet and 
the metallicity of the host: 



P([Fe/H] 



0.03 x 



10 2.o[Fe/H]_ 



(16) 



The relationship is derived from a sample of 850 stars 
with uniform planet detectability in Keck, Lick and 
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Fig. 1. — Upper limits at 95% confidence for each individual clus- 
ter, for both HJ and VHJ period ranges and radii 1.0 Rj and 1.5 
Rj. Solid arrows represent upper limits derived without any rescal- 
ing for the cluster metallicity, and the dashed arrows are rescaled to 
{[Fe/H]) = 0.09 using the FV05 relation between host star metal- 
licity and planet frequency. In cases for which the upper limit is 
formally f p > 1.0, we have plotted the upper limits as f p = 1.0. 



Anglo- Australian Observatory radial velocity surveys, 
and is valid for planet periods of less than 4 yrs, velocity 
amplitudes K > 30 m/s, and -0.5 < [Fe/H] < 0.5. The 
weighted average metallicity of all of the cluster stars 
(not including those from NGC 2362) used in our anal- 
ysis is [Fe/H] =0.09. The cluster NGC 2362 does not 
appear to have any metallicity estimate present in the 
literature, and so we have assigned it the average metal- 
licity of the remainder of our sample. 

Cluster metallicities are subject to large uncertain- 
ties in the literature, and different authors may pub- 
lish metallicities discrepant by up to ~ 0.5 dex 
in [Fe/H]. In the case of NGC 215 8, lite rature es- 
timates range f rom -0.64 ±0-24 (jJanesI I1979D to 
-0.238 ± 0.064 (iTwarog et all fl997h to -0.03 ± 0.14 
Jacobson et alj 120091 ). We adopt the Uacobson et al.l 
2009f l value. Likewise, estimates of the metallicity 



of NGC 6791 typic ally range in the l it erature from 
fFe/ Hl=+0.35-0.45 dGratton et al.1 [2001 lOriglia et al.1 
[2001 lAnthonv-Twarog et al.ll2007t iCarraro et alj |2006) . 
We adopt a metalli c ity fo r NGC 6791 of [Fe/H] = +0.39 
from ICarraro et all (|2006l ). In addition to the significant 
discrepancies among literature values, cluster metallici- 
ties are often quoted without uncertainties. We therefore 
stress that while it is useful to include our knowledge of 
the planet-metallicity correlation our analysis, the result- 
ing upper limits should be considered within the context 
of the uncertainties. 
We derive scaling factors that account for the planet 
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TABLE 3 

Upper limits on the frequency of HJs and VHJs normalized to [Fe/H] = 0.09 











VHJ 






HJ 




Cluster 


S \Fc/R] 




1.0 Rj 


1.2 Rj 


1.5 Rj 


1.0 Rj 


1.2 Rj 


1.5 Rj' 


M37 


0.82 




0.033 


0.031 


0.028 


0.101 


0.094 


0.084 


NGC 188 


0.54 


firm 


0.409 


0.295 


0.123 


3.44 


2.48 


1.03 






marginal 


0.095 


0.077 


0.049 


0.328 


0.265 


0.171 


NGC1245 


0.52 


0.463 


0.327 


0.123 






0.695 


NGC 2158 


0.57 


firm 


0.133 


0.086 


0.016 


0.693 


0.449 


0.082 






marginal 


0.036 


0.026 


0.011 


0.092 


0.066 


0.028 


NGC2362 


1.00 




0.521 


0.369 


0.142 


4.99 


3.30 


0.775 


NGC 6791 


3.93 


firm 


0.021 


0.013 


0.002 


0.088 


0.056 


0.008 






marginal 


0.006 


0.004 


0.002 


0.018 


0.012 


0.005 



metallicity correlation: 

7? ([Fe/H]) 
°[Pe/H] 77(0.09) ' 

which alters the form of Eqn. [10] to 



» < 

p - s 



In (0.05) 



[Fe/H] 



N c iV e 



(17) 



(18) 



since the probability of detecting a planet around a given 
star is enhanced or depreciated by a factor of SVpg/jj] m 

comparison to the average star in the sample. Such a 
scaling means that the null results of metal-rich clusters 
like NGC 6791 become more important in the calculation 
of the combined upper limit. 
Values of -Spg/jj] for each cluster are as follows: 0.82 

for M37, 0.54 for NGC 188, 0.52 for NGC 1245, 0.57 
for NGC 2158, 1.0 for NGC 2362, and 3.93 for NGC 
6791. In the case of NGC 2362, we have assumed that 
it has the average metallicity of the entire sample, and 
f p is therefore unchanged by the inclusion of metallic- 
ity information. The upper limits to the planet fractions 
adjusted for metallicity are given in Table [3] When com- 
bined, these values yield an upper limit on the planet 
frequency of 4.3% for 1.0 Rj HJs, 2.9% for 1.2 Rj HJs, 
0.69% for 1.5 Rj HJs, 1.1% for 1.0 Rj VHJs, 0.72% for 
1.2 Rj VHJs, and 0.16% for 1.5 Rj VHJs, with all values 
quoted here normalized to [Fc/H] = 0.09. 

6. DISCUSSION 

6.1. Comparison to other planet search results 

The frequency of short-period planets has been 
previously e s timat ed in a number of other surveys. 
IGould et al.1 (|200l hereafter G06) used the OGLE III 
microlensing surveys and found the frequency of short- 
period planets in the Galactic field (at 90% confidence) to 
be ~ 0.31% for 3 < P < 5 days and 0.14% for 1 < P < 3 
days ; valid for plan e tary r adii of 1.0 < Rj < 1.25. 
The iCumming etaTl ([20081 hereafter C08) analysis of 
the radial velocity Keck Planet Search results suggest 
that 7/585 of FKG stars have companions with mass 
Msini > 0.1 Mj and periods 3 < P < 5 (see Fig. 5 in 
C08) and 1/585 stars have planets in the p eriod range 
of 1 < P < 3 days. iWeldrake et~aT1 (|2008l ) found up- 
per limits on the frequency at 95% confidence of 0.67% 
for 3 < P < 5 days and 0.096% for 1 < P < 3 days 
in the globular cluster w Centauri. The [Paulson et al.l 



(2004) radial velocity survey of the Hyades open clus- 
ter detected no close in giant planets among 94 stars, 
although upper limits on the planet frequency were not 
calculated. 

In comparing the results of other planet searches to 
the upper limits we calculate for open clusters we must 
be aware of biases inherent in each survey method. Dif- 
ferent survey methods are particularly sensitive to differ- 
ent planetary parameters, and may not always probe the 
same underlying population. One must also be careful to 
account for the effects of metallicity of the population of 
stars surveyed in each study. This is particularly true in 
the case of radial velocity surveys, which tend to be bi- 
ased towards more metal rich stars, and in which planet 
detectability is determined by the mass as opposed to 
the radius of the planet. 

The latter effect is primarily an issue when one wishes 
to compare radial velocity survey results to transit sur- 
veys. The detectability of a planet in a transit survey is 
not a function of its mass so much as its radius, which 
is not the case in radial velocity surveys, in which the 
planet's radius does not pl ay a role in de t ermin ing the 
detectability of the object. [Fortnev et al.l (|2007l ) shows 
that 1.0-1.2 Rj objects can span three decades in mass, 
from nearly stellar to down to super-earth masses, de- 
pending on the composition of the planet. Radial veloc- 
ity surveys, on the other hand, are sensitive to the mass 
of the object, and can span about an order of magnitude 
in radius for different compositions. For most ground 
based transit surveys and typical radial velocity preci- 
sions both methods tend to probe Jupiter-like objects, 
but it is important to note that the populations probed 
in each survey method can differ. For definiteness, we 
will assume that all planets detected in RV surveys with 
Msini > 0.1 Mj have radii > Rj 

Additionally, we must account for the metallicity bi- 
ases inherent to each survey method. To better compare 
our results to those of other surveys we must normal- 
ize our upper limits to the metallicities typical of the 
comparison surveys. FV05 performed a careful anal- 
ysis of stellar metallicities in Keck, Lick, and Anglo- 
Australian Observatory radial velocity surveys, and they 
find that the mean metallicity of 1040 planet search 
stars is 0.09 dex more metal-rich than volume-limited 
sample selected from the same dataset. Furthermore, 
stars that bear planets are 0.226 dex more metal rich 
than the volume-limited sample. The average metallic- 
ity of planet-bearing stars in FV05 is [Fe/H] w 0.142 
which implies a mean metallicity for the RV planet search 
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TABLE 4 
Combined upper limits on f p 



Rj Period w/o [Fe/H] correction ([Fe/H]) = 0.09 a ([Fe/H]} = 0.006 {[Fe/H]) = -0.084 







firm 


marginal 


firm 


marginal 




90% 


1.0 


VHJ 


0.0137 


0.0066 


0.0109 


0.0044 


0.0073 


0.0037 


1.2 


VHJ 


0.0095 


0.0049 


0.0072 


0.0031 


0.0048 


0.0032 


1.5 


VHJ 


0.0031 


0.0024 


0.0016 


0.0013 


0.0011 


0.0006 


1.0 


HJ 


0.0549 


0.0194 


0.0431 


0.0123 


0.0290 


0.0148 


1.2 


HJ 


0.0385 


0.0143 


0.0286 


0.0088 


0.0193 


0.0127 


1.5 


H.I 


0.0139 


0.0068 


0.0069 


0.0036 


0.0046 


0.0024 



a [Fe/H] = 0.09 corresponds to the average metallicity of the cluster stars, [Fe/H] = 0.006 to that of 
stars in RV planet searches, and [Fe/H] = —0.084 to the stars observed in OGLE. 
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Fig. 2. — Combined upper limits utilizing the six selected surveys. The left panel shows upper limits derived without metallicity rescaling 
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stars of [Fe/H] = 0.142 - 0.226 + 0.09 = 0.006 and 
[Fe/H] = 0.142 - 0.226 = -0.084 for the volume-limited 
sample. We therefore adopt [Fe/H] = 0.006 as the aver- 
age metallicity of the stars used in C08 to estimate the 
planet frequency in the Keck survey, which based on a 
target sample similar to that analyzed in FV05. 

In the case of the OGLE III survey, there is only a 
very weak bias towards more metal rich stars relative to 
a volume limited survey (G06). We therefore adopt the 
FV05 average metallicity for the volume-limited sample 
as the mean metallicity of the OGLE stars. Although 
the FV05 volume-limited sample only extends out to 20 
pc while the OGLE survey detects much more distant 
stars, we assume that the FV05 volume-limited average 
metallicity is a fair estimate of the average OGLE metal- 
licity. The uncertainties in the determination of the clus- 
ter metallicities are far larger than those implicit in this 
assumption, and so for our purposes, [Fe/H] = —0.084 
is a sufficiently accurate estimate of the average OGLE 
field star metallicities. 

We cannot better norm alize the Wcldra ke et al I ( 2008) 
and iPaulson et all (|2004fl results. iPaulson et all (2004) 
does not calculate explicit upper limits or detec- 
tion efficiencie s , wh ich makes comparison difficult. 
iWeldrake et al.l ()2008l ) focuses on the globular cluster to 
Cen which has [Fe/H] < -0.5, for which the FV05 rela- 



tion no longer holds. Therefore, although we quote the 
results of these surveys for completeness, we cannot di- 
rectly compare them to our results. 

We can quote our upper limits normalized to the mean 
metallicities of the C08 radial velocity surveys and the 
OGLE transit surveys. For ([Fe/H]) = 0.006, corre- 
sponding to radial velocity surveys, we find upper limits 
of 2.9%, 1.9% and 0.46% for 1.0, 1.2, and 1.5Rj HJs, 
respectively. For VHJs the limits are 0.73%, 0.48%, and 
0.11% for 1.0, 1.2, and l.hRj respectively. For an av- 
erage metallicity or [Fe/H] = —0.084, corresponding to 
the OGLE III surveys, we find upper limits of the planet 
fraction of HJs at 90% confidence (as in OGLE) 1.5% 
of 1.3% and 0.24% for 1.0, 1.2, and 1.5 Rj planets re- 
spectively. Similarly, for VHJs we find limits of 0.37%, 
0.32%, and 0.06% for 1.0, 1.2, and 1.5 Rj planets re- 
spectively (The frequencies quoted in G06 were ~ 0.31% 
and ~ 0.14% for HJ and VHJ respectively for 1.0-1.25 
Rj planets) . Figures [3] and 0] graphically display these 
results. 

For each comparison, our upper limits are consistent 
with the results of other surveys. Our derived upper lim- 
its for 1.0 and 1.2 Rj lie above the short-period planet 
frequencies derived using both the OGLE transit surveys 
and the C08 result for the Keck Planet Search. Another 
way to see this is to ask: how many planets should we 
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expect to have seen in these 6 open cluster surveys given 
the frequency of short-period giants published in G06 and 
C08, given the known detection efficiencies and number 
of cluster stars targeted by each survey? The answer is 
0.86, 1.00, and 5.80 for 1.0, 1.2, and 1.5 Rj VHJs, respec- 
tively, and 0.48, 0.56, and 3.03 for 1.0, 1.2, and 1.5 Rj 
HJs using the G06 frequencies. Using the C08 frequen- 
cies we estimate the surveys should have detected 0.70, 
1.06, and 4.69 planets for 1.0, 1.2, and 1.5 Rj VHJs and 
1.24, 1.86, and 7.75 planets for HJs. The probability of 
a non-detection when one expects to detect 1.86 1.2 Rj 
planets is 15.5%. Although it is possible that these sur- 
veys could yield a non-detection for a planet frequency 
equal to that found in RV surveys, the probablity is rel- 
atively low, suggesting that these transit survey results 
are approaching interesting constraints. However, as we 
discuss in ^6.41 the number of additional observed stars 
needed to significantly improve these existing upper lim- 
its becomes prohibitively large. 

The expected number of planets is only large in the 
case of 1.5 Rj objects, and both the G08 and C06 fre- 
quencies are not directly comparable to planets in this 
radius range. For the OG LE surveys, t his is because 
the upper limits quoted in iGould et al.l (|2006| ) are for 
1.0-1.25 Rj planets; no planets were found with radii 
larger than this, and any limit on 1.5 Rj planets from 
the OGLE survey would be tighter. In the case of the 
C08 result, one must recall that radial velocity surveys 
probe planetary mass, not radius. It is likely the case 
that the planets with Msini > 0.1 Mj detected in RV 
surveys are of smaller radii. The upper limits for a more 
likely average radius of 1.2 Rj for such planets, are not 
in conflict with the C08 frequencies. 

6.2. RV detected planets in open clusters 

We have evidence from other methods of planet 
detection that plane ts do exist in open clusters: 
lLovis fc Mavorl (|2007t ) reported the RV detection of a 
massive planet with a minimum mass of 10.6 Mj on 
a 714 day orbit around a red gian t in t he open clus- 
ter NGC 2423. Likewise. ISato et all (|2007| ) reported the 
discovery of a planet in the Hyades open cluster with 
ui sin i = 7.6 ±0.2 Mj and period 594.9 ±5.3 days, again 
using radial velocities. Although these planets have much 
longer periods than any of those to which these transit 
surveys would be sensitive, they still represent evidence 
that open clusters are not devoid of massive planetary 
bodies. 

6.3. What about all of the other surveys? 

We have chosen a subset of 6 transit surveys from ~ 20 
such projects because these surveys provided us with 
sufficient information to combine their results quantita- 
tively, as described in iJ3J However, this means that we 
have neglected the majority of the photometric surveys 
of open clusters in the literature. Here we attempt to 
ascertain how the addition of these surveys might alter 
our conclusions. 

Although the remaining surveys do not meet our cri- 
teria for selection, we do our best to estimate to what 
degree these additional data affect the upper limits on 
the planet frequency. In most of these surveys the au- 
thors made no attempt to put upper limits on the fre- 
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quency of planets in the fields they observed, often se- 
lected transits by eye, and rarely estimated the num- 
ber of cluster members observed. We estimate from the 
contents of each paper the number of stars with suffi- 
ciently precise photometry for planetary transits to be 
visible, and the number of cluster members likely to be 
within this subset of stars. A list of the additional sur- 
veys and estimates of the number of stars they contribute 
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TABLE 5 
Additional Photometric Surveys 



Cluster 




steirs with 
~ 1% phot. 


cluster stars 


method of membership 
est inictt ion 


Reference 


NGC 2301 


+0.06 


4000 


~ 1500 


estimate from CMD b 


Howell et al. (2005) 


NGC 2632 


+0.14 




~ 150 


WEBDA memberships 


Pepper et al. (2008) 


NGC 2660 


-0.18 


3000 


~ 600 


comparison fields 


von Braun et al. (2004) 
von Braun et al (2005) 


NGC 6208 




5000 


~ 1000 


comparison fields 


Lee et al. (2004) 


NGC 6633 




2000 


~ 600 


by-eye estimate 15 


Hidas et al. (2005) 


NGC 6819 


+0.07 


11500 


~ 700 


colors 


Street et al. (2003) 


NGC 6940 


+0.01 


4400 


~ 100 


most cluster members saturated 


Hood et al. (20051 


NGC 7086 




445 


~ 150 


Besancon model 


Rosvick fc Robb (2006) 


NGC 7789 


-0.08 


2400 


~ 240 


star counts on CCD chips b 


Bramich et al. (2005) 


NGC 6253 c 


+0.36 


? 


? 




Montalto et al. '(2009) 



a Metallicities courtesy of WEBDA: http://www.univie.ac.at/webda/ 

b Even rough estimates of the membership were not mentioned in the the reference, and we estimated membership 
ourselves. In all other cases the references had at least some, however rough, estimate for the number of cluster 
st ars observed. 

c IMont alto et al. (2009) contained no information about the quality of the light curves or details about any 
planetary transit searches. 



to the total sample of cluster stars is presented in Ta- 
ble [5] In total, the 10 additional surveys add ~ 5000 
stars to our sample. Assuming an average detection ef- 
ficiency of 0.5%, and a maximum efficiency of 1.0%, this 
implies that f Pother ~ (5ooo)(o oi) — ^% m the best case 

and fpother ~ (5000)(o.oo5) - 12% on average. The con- 
straints on the population of 1.0 Rj HJ planets from 
the main 6 surveys are the weakest of any considered 
here with f p < 0.055. Using our standard method of 
combining upper limits, the addition of the 5000 stars 
from the other surveys yields a new limit of f p < 0.029, 
which is a ~ 50% tighter constraint. In the case in which 
V t — 0.005, the combined upper limit is f p < 0.037. In 
the VHJ range, the upper limits on f p for 1.0 Rj plan- 
ets decreases by 18% and 10% for detection efficiencies of 
1.0% and 0.5% respectively. At larger planetary radii the 
changes in f p for 1.5 Rj planets due to the addition of 
these 5000 stars is of order ~ 5 — 20%. The conclusions do 
not qualitatively change with the inclusion of metallicity 
information. The mean weighted metallicity of the stars 
in the unused 10 surveys is ([Fe/H]) ~ +0.1 (according to 
the WEBDA database), comparable to ([Fe/H]) w +0.09 
found for the 6 surveys utilized in this paper. 

We expect that by neglecting the observations present 
in these other surveys our derived upper limits may be 
up to a factor of two too high. However, it is highly 
unlikely that all of these surveys have achieved the opti- 
mistic 1.0% detection efficiency, and improbable that all 
have even achieved a 0.5% efficiency. We therefore con- 
clude that neglecting these other surveys does not quali- 
tatively change our primarily conclusion, that the lack of 
detections is consistent with the hypothesis that cluster 
and field stars host the same planet population. We do 
note that, had each of these surveys carefully quantified 
their detection efficiencies, the constraints on the num- 
ber of short-period planets in open clusters could have 
been noticeably tighter. 

6.4. Future transit surveys in open clusters: how many 
more stars are needed? 

It is useful to ask, given our results, how many more 
cluster stars must be observed in surveys with null re- 



sults before the combined upper limit derived using all 
extant surveys becomes inconsistent with the observed 
frequency of planets around field stars. Because of the 
low detection efficiencies and relatively few suitable open 
clusters in our galaxy, it is unlikely that upper limits 
derived from transit surveys in open clusters will be in- 
consistent with the results of C08 and G06 in the near 
future, even if it is the case that the frequency of plan- 
ets in open clusters is significantly lower than that of 
the field. For the C08 result, a combined upper limit (at 
95% confidence) would be inconsistent with the 68% con- 
fidence lower bound on the C08 frequency of the fraction 
of stars with planets was found to be f p < 0.0076 for 
1.2 Rj planets in the 3 < P < 5 day period range. The 
six-survey upper limit is f p < 0.019, which implies that 
~ 48, 000 more cluster stars must be observed (with null 
results) at a 0.5% detection efficiency before the limit 
is in conflict with C08 (currently, ~ 13000 stars have 
been observed at ~ 0.5% efficiency). This number drops 
to 9000 with a detection efficiency of 1.0%. Given that 
there are about 1200 known open clusters in the galaxy 
(WEBDA), only a handful of which are rich enough to 
contain thousands of stars, the only way to make the 
upper limits on the planet frequency in open clusters de- 
rived from transit surveys competitive with field surveys 
is to drastically increase the detection efficiency of the 
surveys and number of clusters observed. When we ask 
how many stars would need to be observed to be com- 
petitive with the G06 field star planet frequency, the re- 
quirements are even more strict. Surveys would need to 
cover ~ 387, 000 stars at 0.5% detection efficiency before 
the upper limits, at 90% confidence, would be inconsis- 
tent with the 68% confidence lower bound on the G06 
frequency. When we consider the upper limits on VHJs 
even more additional stars are required. 

7. CONCLUSION 

We conclude that the current upper limits of the HJ 
and VHJ frequencies determined using the null results 
of transit surveys in open clusters do not suggest a sig- 
nificant difference between the frequency of planets in 
open clusters and the field. Open clusters remain, to the 
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best of our knowledge, a viable and useful target for ex- 
oplanet transit surveys, and we do not yet have enough 
data to discern whether the environments of open clus- 
ters have any noticeable effect on planet formation and 
survival. We recommend that any future surveys care- 
fully quantify any null results, since the combination of 
many such outcomes has the potential to better constrain 
the frequency of short-period planets and answer inter- 
esting questions about the formation of planets in stellar 
clusters. 
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